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' Context. Thanks to the improved angular resolution of modern telescopes and kinematic models, the existence of 

\Q , supermassive black holes (SMBHs) in the inner part of galaxies, regardless their morphology and nuclear activity, has 

been established on quite solid grounds. A possible correlation between the mass of SMBHs (Af.) and the evolutionary 
state of their host galaxies is expected and is currently under a heated debate. 

Aims. Based on the recent 2D decomposition of 3.6 fim Spiter /IRAC images of local late- and early-type galaxies with 
M, measurements, we investigated various scaling laws, studying what the best predictor of the mass of the central 

O black holes is. In particular, we focused on the M. — Mqo 1 law, that is the relation between the mass of SMBHs and 
the kinetic energy of random motions of the corresponding host galaxies, Mq is the mass and a the velocity dispersion 
(— I ' of the host galaxy (bulge). 

Methods. In order to find the best fit for each of the scaling laws examined, we performed a least-squares regression 
of M, on x for the considered sample of galaxies, x being a whatever known parameter of the galaxy bulge. For this 
O ' purpose, we made use of both the linear regression LINMIX_ERR and FITEXY methods. 

Results. Our analysis shows that M, — Mg<t 2 law fits the examined experimental data successfully as much as the other 
known scaling laws (all correlations have similar intrinsic scatters within the errors) and shows a value of x 2 (estimated 
by FITEXY) better than the others, a result which is consistent with previous determinations at shorter wavelengths. 
This means that a combination of a and Mg (or R e ) could be necessary to drive the correlations between M, and other 
bulge properties. This issue has been investigated by a careful, although not fully conclusive, analysis of the residuals 
of the various relations. 

Conclusions. In order to avoid rushed conclusions on galaxy activity and evolution, the indirect inferring of the masses 
of the supermassive black holes from the kinetic energy of random motions via the M, — Mgct 2 relation should be 
considered, especially when applied to higher redshift galaxies (z > 0.01). This statement is suggested by a reanalysis 
of the Sloan Digital Sky Survey (SDSS) data used to study the black hole growth in the nearby Universe. By adopting 
the M, — Mgct 2 relation instead of the M. — a relation, a radio-quiet/radio-loud dichotomy appears in the SMBH 
mass distribution of the corresponding SDSS early-type AGN galaxies. 

Key words, black hole physics - galaxies: general - galaxies: kinematics and dynamics - galaxies: statistics - galaxies: 
active - galaxies: evolution 



K> , 1. Introduction well as the structural para meters of the host galax- 

j_j ■ ies (bulges) (|Ma,gorrian et al.l IT998I: iTremaine et aD 120021: 

& j Scaling laws in galaxies are important to describe the mech- [ Marconi fc Huntll2003l: iGebhardt et al.ll2003t iHaring fc Rbd 

anisms for the initial formation of the first galaxies as 2004 1: lAller fc Richstond 120071: iGrahaml 12008 alibi: lHull200S 



well as for their cosmic evolution. They also apply quite Kisaka et al.l 2008: Fc oli fc Mancinil 20091: iGiiltekin et aD 

well to high redshifts, implying that the interaction and l2009al h 

merging process es must, at least on average, preserve them Thanks to these catalogues, the astrophysical com- 

(|Schneider||2006|) . munity identified a large number of scaling laws, 

One of the recent successes in extragalactic astronomy in which the mass of SMBHs correlates with sev- 

was the discovery that both early- and late-type galax- eral properties of the host spheroidal component], 

ies, close (< 100 Mpc) to the Milky Way, host a super- such as for instance bulge luminosity, mass, effective 

mas sive black hole (SMBH; M, > 1O 6 M0) at their cen - radius, central potential, dynamical mass, concentra- 

ter (|Kormendv fc Richstond 119951: iRichstone et~aD Il998l) . 

Subsequently, based on more and more precise experi- i Here we use the terms bulge or spheroid to mean the 
mental data, it was possible to draw up many galaxy spheroidal component of a spiral/lenticular galaxy or a full el- 
data sets, containing measures of the SMBH masses as liptical galaxy. 
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f 995; Magorrian ct al 
Gephardt et al.ll2000l 

Wandelll2002t iTremaine et al.ll2002UMarconi fc Huntll200_ 
Hari ng fc Rixl l2004 iGraham fe Driverl 120051: iFeoli fc Mel 



1998; iFerrarese fc Merrittl 2000; 



Laorll2001fclMerritt fc Ferrareseil200ll 



tion, Sersic index, bin d ing e n ergy, etc. (iRichstone et all bulge properties: kin etic energy M. — M+a 2 , momentum 

parameter M.-M+o (ISoker fc Meironll201ll) . and the M. - 
R e <7 3 law. In order to have a comprehensive stu dy, we also 
reanal yzed the relations already investigated by ISani et all 
([201 If ). This allows us to have a rapid comparison among 
the various relationships and to find what is the tightest. 

Secondly, we examined if a simple one-to-one correla- 
tion b etween, e.g . . M, and a is an exhaustive description 
of the ISani et all ( 20111 ) data, or if we have to consider an 
additional dependence on a second parameter such as R c or 
M+. In order to study the ex istence of such a depend ence, 
we used both the approach of IMarconi fc Hunt] (|2003f ). who 
investigated the correlation of the r esidual of the M« — a re - 
lation with R e and M*, and that of lHopkins et al.l (|2007aD . 
who considered the correlations between residuals at fixed 



2005t lAller fc Richstond 120071: IHopkins et al.l l2007bl) . 
Another s caling law has been rec ently proposed, the M, — 
R e a 3 law (|FeoU fc Mancinill2011l) . which is based on a pure 
theoretic al framework, and is on the wake of the numerical 
results of IHopkins et al.l (|2007al) . 

Relation s between M, and X -ray luminosity, radio- 
lumi nositv dGultekin et all l2009bl) . momentum paramc- 



2011), number of globul ar clusters 
20101: ISnvder et all 1201 ll ) have also 



ter (ISoker fc M eiron 
(jBurkert fc Tremaind 

been presented, whereas the correlation wi th the dark mat- 
ter halo is still a matter of controversy dFerraresd [2002; 
Baes et al.ll2003HKormendv fc BendeHl201lUGrahamll201l[ 



Volonteri et alJl20lH iBellovarv et al.ll2011|) . 



All these scaling laws have led to the belief that SMBH 
growth and bulge formation regulate each other ( Ho 2004) , 
even if it is difficult to understand the fundamental na- 
ture of_tlie_TOr£eiationsl2etween SMBHs and host proper- 
ties ([Jahnke fc M actio 20111) , also because all such relations 
depend critically on the accuracy of the publis hed error es- 
timates in all quantitie s under consideration ([Novak et all 
120061: lLauer et alll2007f ). 

Just as the Faber-Jackson relation, the "traditional" re- 
lation between the SMBH mass and bulge velocity disper- 
sion a or stellar mass M* should be projections of the same 
fundamental plane relation. Current observations require a 
correlation of the form M, cx M"a^ over a simple correla- 
tion w i th either a or M+ at > 3a confiden ce ([Hopkins et all 
l2rjTJ7llHopkn^l2rMlMarulli et al.|[2008h . 

Actually, another competitive correlation, as opposed 
to the popular M. — M* and M. — a, between M. and 
the kinetic energy of random motions of th e correspond- 
ing bulges, i.e. M*o 2 has been advanced ([Feoli fc Meld 
120051: iFeoli fc Mancinill2009D . This relation has also a plau- 
sible physical interpretation that resembles the H-R dia- 
gram: the mass of the central SMBH, just like entropy, 
can only increase with time or at most remain the same 
but never decrease; M. is therefore related to the age of 
the galaxy. On the other hand, the kinetic energy of the 
stellar bulges directly determines the temperature of the 
galactic system. The goodness of the M, — M+a 2 relation, 
as a predictor of the SMBH mass in the center of galax- 
ies, has been already tested, with clear positive results, 
over three independent galaxy samples and in the frame- 
work of the ACDM cosmology, using two galaxy forma- 
tion models based on the Millennium Simulation, one by 
Bower et al.l (120061) ( the D urham model) and t he other by 



De Lucia fc Blaizotl ([2007D (the MPA model) (|Feoli et all 
20111) . 

Recently, ISani et all ([201 ll) have presented a mid- 
infrared investigation of the scaling relations between 
SMBH masses and some of structural parameters (luminos- 
ity, mass, effective radius, velocity dispersion) of the host 
spheroids in local galaxies, based on Spitzer /IRAC 3.6 /im 
images of 57 galaxies of different morphological types with 
M, measurements. Their results were consistent with the 
above mentioned determinations at shorter wavelengths. 

The prese nt work has severa l aims. First of all, we used 
the data set of lSani et al.l (|2011f ) and completed their study, 
by analyzing the other known scaling relations that involve 



Finally, by considering the same Sloan Digital S ky 
Survey (SDSS) data set, used bv lSchawinski et aTT(|2010l ) to 
study the role of SMBH growth in the evolution of normal 
and active galaxies, we discuss the possible consequences of 
the use of the scaling law M. — M±a 2 in the place of M. — a 
in inferring the SMBH masses. 

The paper is structured as follows. In ^Hwe describe the 
galaxy sample examined and the fitting procedures per- 
formed to find the best-fitting lines for each of the rela- 
tionships considered. In ^3]we reported the fitting param- 
eters and showed the plots of the various scaling relations. 
The issue related the existence of a SMBH fundamental 
plane is tackled in 2) In order to emphasize how the con- 
clusions could be different using a scaling law rather than 
another one, in fl 5] we remake the analysis performed by 
ISchawinski et al. (2010) on SDSS data. After that we draw 
our conclusions in |J6J 



2. Data analysis 

An important fact that emerg e s from the Spitzer /IRAC 
images, analyzed bv ISani et all (|201ll ). is that the 3.6 /Ltm 
luminosity proved to be a very good tracer of stellar mass. 
Actually, this is a crucial information in order to perform a 
proper 2D photometric decomposition of the galaxy compo- 
nents (disk, bar, bulge, etc.), which is useful to investigate 
the interplay between central SMBHs and the evolutionary 
state, luminosity and dynamics of their host galaxies. 

Thanks to this bulge-disk decomposition, the authors 
were also able to identify in their sample 9 disk galaxies 
that host a pseudobulge0, i.e. Circinus, IC2560, NGC1068, 
NGC3079, NGC3368, NGC3489, NGC3998, NGC4258, and 
NGC4594. Constructing the correla tion between M« and 
the bulge 3.6 ^m luminosity Li fi v,„i. ["Sani et all (|201ll) no- 
ticed that four of these nine are consistent with classical- 
bulge galaxies, whereas the other 5 are outliers at more 
than 4cr below their linear regression. These galaxies are 
Circinus, IC2560, NGC1068, NGC3079, and NGC3368. 
This result is consistent with the fact that the M. — a 
relation for pseudobulges is different from the relation i n 
the classical bulges at a significance level > 3cr |Hu 2008), 
and that at a fixed bulge mass, M. in pseudobulges are on 
average more than one mag nitude smaller than the ones 
in classical bulges ([Hul I2009T ) . Moreover the elliptical-only 



2 Essentially, a pseudobulge is a bulge that shows photomet- 
ric a nd kinematic evidence for disk-like dynamics llKormendvl 



1993). 
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galaxies, and the non-barred galaxies, define tighter rela- 
tions with less scatter and a reduced slo pe than the one 
obtained when usin g a full galaxy sample (|Grahamll2008bt 
iGraham etafl 1201 lh. For these reasons, we neglected from 
the sample of ISani et al.l (|2011h the above mentioned 5 
galaxies and thus considered a more consistent sample of 
N = 52 galaxies, which is therefore formed by 24 ellipticals, 

3 dwarf ellipticals (dEs), 11 lenticulars, 6 barred lenticulars, 

4 spirals, an d 4 barred s pirals . On the contrary the fits in 
the paper of ISani et al.l (|201lh considered only 48 objects 
excluding all the nine pseudobulges. 

By using this sample (all t he pa rameters are reported on 
Table 2 and 3 of ISani et al.l (|201lD ). we investigated what 
the relationship that best predicts the black hole mass is. 
In particular, the relations that we want to study can be 
written in the following form 



log 10 M, = b + m log 10 x, 



(1) 



where m is the slope, b is the normalization, and x is a 
parameter of the host bulge. Equation ([T]) can be used to 
predict the values of M, in other galaxies once we know the 
value of x. In order to minimize the scatter in the quan- 
tity to be predicted, we have to perform an ordinary least- 
squares regression of M, on x for the considered galaxies, of 
which we already know both the quantities. We considered 
error bars in both variables and, to simplify the analysis, 
we make all of them symmetric about the preferred value 
by averaging the size of the upper and lower 1 a error bars 
so that xZ7 becomes x±(h+l)/2. To obtain the parameters 
of the fits (m and b), we adopt the following three different 
fitting methods. 



1) The linear regression routine FITEXY ()Press et al.l 
1992) for the relation y = b + mx, by minimizing the x 2 



N 

E 



mxi 



^ (Ay,) 2 +^ 2 (A^) 2 ' 



(2) 



The most efficient and unbiased estimate of the slope 
is obtained when the fitting method incorporates the 
residual variance, also known as intrinsic scatter Sq, 
which is that part of the variance which can not be at- 
tributed to specific causes (jNovak et alJ l2~006). So, if the 
reduced x? = X 2 /(N — 2) of the fit is not equal to 1, 
we normalize including the suitable value of Sq in the 
Eq.© to obtain 



1 N 



(yi-b- mxif 



N - 2 



-f (Ay,) 2 + e 2 + m 2 (Axi) 2 



(3) 



Lastly, it is possible to obtain an estimate of the la 
error bar on £o by adjusting it until the x 2 is equal to 
l + (2/A^) 1 /2. 

2) The linear regr e ssion FITEXY method as modified by 
iTremaine et a l. (2002), where the measurement errors 
of the dependent variable and the intrinsic scatter are 
added in quadrature, adjusting e and refitting until the 
reduced x 2 °f the fit is equal to 1. 

3) The Bayesia n linear regression routine LINMIX_ERR 
(|Kellvll2007f ) to determine the slope, the normalization, 
and the intrinsic scatter of the relationship 



This routine approximates the distribution of the in- 
dependent variable as a mixture of Gaussians, bypass- 
ing the assumption of a uniform prior distribution on 
the independent variable, which is used in the deriva- 
tion of x 2 -FITEXY minimization routine. Since a direct 
computation of the posterior distribution is too compu- 
tationally intensive, random draws from the posterior 
distribution ar e obtained u sing a Markov Chain Monte 
Carlo method (|Kellvll2007D . 

In order to b e consistent with the results reported by 
ISani et al.l (|201ll) . we use the LINMIX.ERR as the favorite 
routine to calculate the fitting parameters of the SMBH- 
bulge scaling relations. 

In Eq. ((TJ in place of x we considered the following quan- 
tities: a, Af dy „, M*, M dyn a, M+a, M dv n a 2 , M+a 2 , R e a 3 



M(jy n is the bulge dynamical mass, which ISani et al.l (| 2 1 lh 
assumed dominated by stellar matter with a negligible con- 
tribution of dark matter and gas, and computed as: 



M dyn = kR e a 2 /G, 



(5) 



log 10 M.=b+ m log 10 (a;) + e . 



(4) 



where G is the gravitational constant, while the factor k 
was fixed by the authors to be equal to 5, in agreement 
with ICappellari et al.l (|2006|) . On the contrary, the stellar 
mass M t of each galaxy has been estimated by combin- 
ing the bulge 3.6 /im luminosity with the galaxy mass-to- 
light ratio (M/L); the values are reported in Tablc[T] (Sani 
fc Ma rconi, private communication). We follow ISani et al.l 
()201lD in this choice, even if we would have preferred an 
estimate of masses by means of the Jeans equatio n or 
Schwarzschild methods as in (jFeoli fc Mancinill2009f ). be- 
cause the use of Eq. (O makes the the relations M. — M dyn a 
and M, — R e a 3 , which are deduced from different theoret- 
ical contexts, practically equivalent. 

3. Results 

Referring to the sample of 52 galaxies, in Table [5] we col- 
lected the parameters of the fits obtained thanks to the 
three above mentioned fitting methods for the various re- 
lations that we analyzed, together with the corresponding 
values of the x" '■> the intrinsic scatter £o, and the Pearson 
linear correlation coefficient r. In Figures 1 — 4, we reported 
the relations in log-log plots (we associated a particular 
marker to each galaxy according to its morphological type). 
The best-fitting lines are also shown for each diagram. 

The relations between M m and the corresponding kinetic 
energy, momentum parameter, velocity dispersion, galaxy 
mass and i? e f 3 , fitted with a Bayesian approach to linear 
regression are: 

log 10 M. = (5.30 ± 0.26) + (0.63 ± 0.05) log 10 [M dyn <r 2 /M c 2 ] ; 

(eo =0.30 ±0.16) 

log 10 Al, = (5.18 ± 0.27) + (0.65 ± 0.06) X log 10 [Af* cr 2 /M & c 2 ] ; 

(eo = 0.30 ±0.16) 

log 10 M. = (8.29 ± 0.05) + (3.95 ± 0.31) X log 10 [cr/200kms- 1 ] ; 

(e = 0.31 ±0.16), 

log 10 M. = (3.15 ± 0.47) + (0.71 ± 0.07) X log 10 [Re a 3 /cGM Q ] ; 

(eo = 0.33 ±0.17), 
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Table 1. The values of the stellar mass and la error for 
each galaxy of the sample (Sani & Marconi 2011, private 
communication.) 



Galaxy 






-) 


Circinus 


10.26 


(0.020 


0.023) 


IC1459 


12.06 


(0.130 


0.160) 


IC2560 


10.94 


(0.260 


0.570) 


IC4296 


12.13 


(0.190 


0.320) 


NGC221 


8.814 


(0.230 


0.450) 


NGC524 


11.74 


(0.080 


0.095) 


NGC821 


11.50 


(0.100 


0.140) 


NGC1023 


10.62 


(0.130 


0.160) 


NGC1068 


11.33 


(0.023 


0.023) 


NGC1300 


10.37 


(0.140 


0.190) 


NGC1316 


11.96 


(0.047 


0.047) 


NGC2549 


10.07 


(0.059 


0.071) 


NGC2748 


10.18 


(0.035 


0.047) 


NGC2778 


9.766 


(0.071 


0.071) 


NGC2787 


9.980 


(0.210 


0.380) 


NGC2974 


11.14 


(0.059 


0.071) 


NGC3031 


11.00 


(0.130 


0.170) 


NGC3079 


11.02 


(0.035 


0.035) 


NGC3115 


10.99 


(0.047 


0.047) 


NGC3227 


10.90 


(0.190 


0.320) 


NGC3245 


10.49 


(0.083 


0.100) 


NGC3368 


10.74 


(0.047 


0.047) 


NGC3377 


10.63 


(0.083 


0.095) 


NGC3379 


11.20 


(0.071 


0.083) 


NGC3384 


10.12 


(0.035 


0.035) 


NGC3414 


11.02 


(0.190 


0.320) 


NGC3489 


10.07 


(0.170 


0.410) 


NGC3585 


11.17 


(0.230 


0.450) 


NGC3607 


11.46 


(0.100 


0.130) 


NGC3608 


11.26 


(0.095 


0.110) 


NGC3998 


10.36 


(0.023 


0.023) 


NGC4026 


10.43 


(0.083 


0.095) 


NGC4151 


10.67 


(0.160 


0.240) 


NGC4258 


10.98 


(0.035 


0.035) 


NGC4261 


11.63 


(0.071 


0.071) 


NGC4374 


11.87 


(0.023 


0.023) 


NGC4459 


10.70 


(0.047 


0.047) 


NGC4473 


11.87 


(0.023 


0.023) 


NGC4486 


11.94 


(0.023 


0.023) 


NGC4486A 


10.06 


(0.035 


0.035) 


NGC4552 


11.13 


(0.035 


0.035) 


NGC4564 


10.68 


(0.059 


0.071) 


NGC4594 


11.23 


(0.059 


0.071) 


NGC4596 


10.81 


(0.023 


0.023) 


NGC4621 


11.49 


(0.047 


0.047) 


NGC4649 


11.69 


(0.059 


0.071) 


NGC4697 


11.33 


(0.035 


0.047) 


NGC5077 


11.75 


(0.100 


0.140) 


CenA 


11.31 


(0.110 


0.150) 


NGC5576 


11.39 


(0.071 


0.071) 


NGC5813 


11.93 


(0.083 


0.095) 


NGC5845 


10.64 


(0.011 


0.011) 


NGC5846 


11.57 


(0.023 


0.023) 


NGC6251 


12.36 


(0.083 


0.095) 


NGC7052 


12.06 


(0.071 


0.083) 


NGC7457 


9.635 


(0.300 


0.910) 


NGC7582 


11.15 


(0.170 


0.270) 



log 10 M m = (2.64 ± 0.53) + (0.71 ± 0.07) x log 10 [M dyn a/M Q c] ; 

(eo = 0.33 ±0.17), 

log 10 M. = (2.47 ± 0.54) + (0.73 ± 0.07) X log 10 [M* er/M c] ; 

(eo = 0.34 ±0.17), 



log 10 Mm = (8.18 ± 0.06) + (0.80 ± 0.09) x (log 10 [Af dyn /M ] - ll) ; 

(eo = 0.38 ± 0.19), 

log 10 M m = (8.15 ± 0.06) + (0.80 ± 0.09) x (log 10 [Af*/Af Q ] - 11) ; 

(e = 0.40 ±0.19). 

By inspection of Table [3J it is possible to note that all 
the fitting parameters are consistent with previous deter- 
minations from the literature at shorter wavelengths. The 
3.6 /im Mm — MdynO -2 relation looks to be slightly preferable 
compared with the other M, — bulge relations, especially if 
one considers the values of \r estimated by FITEXY (col- 
umn 5). However, since all correlations have similar intrin- 
sic scatters within the errors (column 6), we cannot conclu- 
sively determine what the best one is. 

4. A possible fundamental plane for supermasssive 
black holes 

By analyzing a sample of 27 galaxies, which are deemed 
to have "secure" SMB H and bulge mass measurements, 
iMarconi fe Hunt! (|2003l) were the first to note that M, is 
significantly correlated both with a and with R e . Plotting 
the residuals of the M, — a correlation against R c , they 
concluded that a combination of a and R c was nec- 
essary to drive the correlations between M, and other 
bulge prop erties. This top i c was then theoretically inves- 
tigated by iHopkins et al.l (|2007al ) by simulations of ma- 
jor galaxy mergers, which defined a fundamental plane 
(FP), analogous to the FP of elliptical galaxies, of the 

form A!, oc i? 1/2 cr 3 or M. oc Ml /2 g 2 , whe r e M+ 
is the bulge stellar mass, and by iMarulli et al.l ([2008) 
who found M, o c (M ^cr 2 ) 07 !^! Moreover, the sample of 
M arconi fc Hunt! (|2003j ) was reanalyzed by IHopkins et al.l 
(2007a), who found that the observations define a FP that 
should be preferred over a simple relation between SMBH 
and any of a, Md yn , Af*, or R e alone at > 3cr (99.9%) 
significance. 

However, lAller fc Richstond (|2007l) noticed that the 
Mm — M,(a) residuals for their sample o f 23 galaxies did 
not in dicate the combination suggested bv lMarconT fc Huntl 
(2003)- The evidence of a correlation between the residuals 
and the effective radius i s obtained by co nsidering only spi- 
ral and lenticular bulges. [Graham! ((2008b) reached the same 
result studying a sample of 40 galaxies. In particular he 
found that the barred galaxies are responsible for much of 
the trend between the M. — a residuals and R c , whereas the 
elliptical galaxies alone do not provide substantial support 
for the existence of a FP plane for SMBHs. The analysis of 
ISani et ahl (|2011l) does not confirm the existence of a FP. In 
fact, comparing the residuals of M, — a for bulges with their 
effective radius, they did not found any significant correla- 
tion either for the entire sample (r = 0.29), or excluding 
barred galaxies and/or pseudobulges (r = 0.20 — 0.29). 

Plotting the dependence of the residua l of the M, — 
a relation on R c and M*, the results of IHopkins et al.l 

Another e f fort in this sense has been performed by 
iGiiltekin et al.l l|2009bl ). who analyzed the relationship among 
X-ray luminosity, radio luminosity, and mass of a sample of 
SMBHs, identifying a FP that can be turned into an effective 
SMBH mass predictor. 
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Table 2. Regression results for logAf. = b ± to log a; with a sample composed by 52 galaxies. Column (1): scaling 
relations. Column (2): linear regression methods. Columns (3)- (4): the regression coefficients, the intercept b and the 
slope to. Column (5): the reduced \ 2 computed by FITEXY. Column (6): the intrinsic scatter. Column (7): the Pearson 
linear coefficient. 



Relation 


Method 


b 


m 


xi 


£0 


r 


M. - 


M dyn a 2 


LINMIX_ERR 

FITEXY 
FITEXY.T02 


5.30 ±0.26 
5.03 ±0.12 

5.31 ±0.24 


0.63 ±0.05 
0.69 ± 0.03 
0.63 ±0.05 


3.85 


0.30 ±0.16 
0.30 ± 0.04 
0.29 


0.85 


Mm ~ 


- M+o 2 


LINMIXJ3RR 

FITEXY 
FITEXY.T02 
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Fig. 1. Mm — MdvnC 2 ( left) and M. — M+a 2 (right) relations for the sample of N — 52 galaxies extracted from the data 
set of [S ani et all (|2011[ ). The symbols represent elliptical galaxies (red ellipses), lenticular galaxies (green circles), barred 
lenticular galaxies (dark green circles), spiral galaxies (blue spirals), barred spiral galaxies (dark blue barred spirals), and 
dwarf elliptical galaxies (orange ellipses). The black lines are the lines of best fit. 



(20071: ) are clearly in conflict with that of the other au- 
thors (Aller fc Richstond 120071: iGrahaml [20085 ISani et al.l 
[20rl . Of course, the explanation of this difference have to 
be found in the differ ent analysis approach. In particular, 
Hop kins et al.l (|2007bl) considered the correlations between 
residuals at fixed u, and not simply the correlation between 
the residual of M. — a and the actual value of R e or M*. 



As stressed bv lHopkins et al.l (|2007b[ ). if we were to do the 
latter, we miss the significance of any real residuals: "the 
slope recovered (i.e., the inferred dependence of M t on R e ) 
is severely biased towards being too shallow for any nonzero 
dependence on R c , and in only 1% of cases will such a 
method recover a slope similar to the true intrinsic corre- 
lation. Looking at the significance of the residuals in this 
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F ig. 2. M, — Mdvn (left) and M, — M* (right) relations for the sample of N = 52 galaxies extracted from the data set 
of lSani et all (|201lh . The symbols are the same as in Figure [T] 
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F ig. 3. M, — a2nn ( left) and M. — i? c 3 (right) relations for the sample of N — 52 galaxies extracted from the data set 
of lSani et al.l (|201 lh - The symbols are the same as in Figure [1] 



space, it is clear that this projection biases against detect- 
ing any significant residual dependence on R c or M*" . 

Here we used both approaches performing the analy- 
sis on both the 57 and 52 galaxy samples. The results are 
reported in Fig. [5] and Fig. [SJ where we highlighted the po- 
sition of pseudo-bulges, together with the values of the in- 
trinsic scatter and the Pearson linear coefficient. The best- 
fitting lines have been obtained through the LINMIX_ERR 
routine. In particular, Fig . [5]sho ws the results obtained us- 
ing the iMarconi fc Hunt! (|2003|) ap proach, whereas Fig, p 
shows the results obtained using the Hopkins et al.l (l2007bfT 
appro ach; Fig. [5^,b and |6^i,b refer to the full ISani et all 
(2011) sample, whereas Fig.[SJ:,d and|BJ:,d refer to the most 
consistent sample of 52 galaxies. 



Both the approaches returned similar results, but it is 
interesting to note how the values of e Q and r slightly im- 
prove moving from the full sample to that of 52 galaxies. 
This fact suggests that the choice of the galaxy sample is 
critical in order to get reasonable results. However, the re- 
sults reported in Table [5] and the above analysis of the 
residuals, applied to th e 52 galaxy sample, are not decisive 
to confirm the result of Hopkin s" et al.l (|2007bh . that is the 
M. - M*ct 2 is preferred over a simple relation between M m 
and any of a or M„ alone. 
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5. Inferring the mass of black holes indirectly in 
high-redshift galaxies 

Scaling relations between astrophysical quantities always 
hide fundamental driving mechanisms. An important step, 
in order to understand these physical mechanisms, is to 
identify where the scaling laws apply and their nature. 
Without this information it is hard to say what the best 
scaling law, linking the mass of the SMBHs with the right 
parameter of the hosting bulges, is. For instance, if we use 
a whatever scaling law in order to infer the masses of the 
SMBHs located in the center of high-redshift galaxies, and 
then we use them to study galaxy-evolution trends, we 
could draw incorrect or misleading state ments. From this 
point of view, the case of the paper of ISchawinski et all 
(|2010T) is emblematic. These authors used SDSS data and 
visual classification of morphology from the Galaxy Zoo 
projectQ to study black hole growth in the nearby Universe. 
They selected all galaxies with SPSS spectra classified 
as GALAXY ([Strauss et al.ll200lD in the redshift interval 
0.02 < z < 0.05; from this parent sample of 47675 they se- 
lected a small (~ 2%) sub sample of 942 narrow-line Active 
Galactic Nuclei (AGN), excluding broad-line AGN, highly 
obscured AGN, and LINERs. Then they inferred the masses 
of the SMBHs indirectly from the stellar velocity disper- 
sion at the effective radius via the M, — a relation usin g 
the slope and the normalization of iTremaine et ail ([2002) , 
(b = 8.13, m = 4.02). Finally, they reported the distribu- 
tion of inferred SMBH masses, plotting only objects where 
the measured a was greater than 40 km s - 1 , corresponding 
to log 10 (Af.) - 5.3 ([Tremaine et al.ll2002ft . 

Using the s ame galaxy catalogue (Table 3 of 
ISchawinski et all (|2010f ) , and a private communication from 
Schawinski, 2011), we inferred the masses of the corre- 
sponding SMBHs indirectly both from the velocity disper- 
sion (via the M. — a relation) and from the kinetic energy of 
random motions (via the M, — Md yn cr 2 relation), using the 



4 www.galaxyzoo.org 



slopes and normalizations obtained by the LINMIX_ERR 
routine and reported in Table [21 In Fig. [7] and Fig. [5] we 
plotted the distributions of inferred SMBH masses from 
both the AGN (colored) and the normal galaxy (white) 
populations, both for the entire sample and split by host 
morphology. These histograms clearly highlight the differ- 
ent distributions of black hole growth in agreement with 
the parent population and with the adopted scaling law. 

Co mparing our Fig. [7] with Fig. 8 of ISchawinski et all 
(|2010l ). we find the same picture: AGN early-type galax- 
ies have lower black hole masses than normal early-type 
galaxies, whereas for late-type galaxies there is an oppo- 
site trend. In early-type galaxies, it is preferentially the 
galaxies with the least massive black holes that are ac- 
tive; instead, in late-type galaxies, it is preferentially the 
most massive black holes that are active. The median val- 
ues of the SMBH masses of both early- and late-type 
AGN host galaxies are still very similar to each other 
(8.38 x 10 6 M Q and 1.17 x 1O 7 M ), but higher with respect 
to the values (2.81 x 1O 6 M and 4.27 x 10 6 A/ Q ) reported 
bv ISchawinski et al.l (|2010l ) due to the different slope used 
in the two analysis. 

On the other hand, if we compa re our Fig. [8]with Fig. [Jj 
or Fig. 8 of Schawi nski et al.l ([2010D . it is possible to observe 
the same trend of the SMBH activity in early- and late- 
type galaxies, but the distributions are much more crowded 
and peaked. As a matter of fact, we found all the distri- 
bution peaks and cut off at greater values of the SMBH 
mass (the median values of the SMBH masses of early- 
and late-type AGN host galaxies are 2.85 x 10 7 Mq and 
3.95 x 1O 7 M ) respectively. 

Another different characteristic, which is quite evident, 
is that the early-type AGN galaxies have a two-peak dis- 
tribution (see upper-right panel in Fig. [5]). This feature is 
still slightly present in the indeterminate-type AGN (Fig. 
Et), whereas it disappears in the late-type AGN (Fig. jSJi). 
It looks like there were two different kinds of early-type 
AGN: one equipped with black holes of low mass (around 
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Log 10 R,(kpc)] Log 10 [M,/M o J 

Fig. 5. Residuals of M. — a as a function of the host galaxy effective radius R e (left) and stellar mass M* (right), for 
the full sample of 57 galaxies (a-b), and for a more consistent sample (see text) of 52 galaxies (c-d), respectively. The 
corresponding intrinsic scatter and Pearson linear coefficient are reported in the upper-left corner of each plot. The 
symbols are the same as in Figure [TJ Galaxies inside the black circles are pseudo-bulges. 



3.2 x 10 6 Mq), and the other one with black holes of larger- 
mass (around 2.5 x 1O 7 M ). 

Actually, a similar double-pe ak figure has been found 
bv lCapetti fc Balmave rdc (2006) plotting the distributions 
of radio-quiet and radio-loud AGN galaxies selected from 
the HST and Chandra archival data. So the double peak in 
the early-type AGN (Fig. \E}p) could be easily explained by 
the radio-quiet/radio- loud dichotomy. 

As already noted bv lChiaberge et a l. (2005), all radio- 
loud AGN are associated with SMBH masses > 1O 8 M , 
whereas most of the radio-quiet population has lower 
SMBH masses. A simil ar res ult has been recently found b y 
iBaldi fc Capettll (l2010h and IChiaberge fc Marconil (1201 111 , 
and is in agreement with what we found in Fig. |BJd. Since 
radio-galaxie s are almost universally found hosted by ellip- 
tical galaxies (jUrrv fc Padovanill995F) , this explains why we 



did not see any double peak in the late-type AGN galaxies 

(Fig. EH). 



It would be interesting to examine a larger sample of 
AGN in order to understand if this double peak is real or 
caused by a too small sampling, but this is beyond the scope 
of this paper. What we want to show here is that if we use 
the M, — Mdyncr 2 instead of the M, — a relation, we obtain 
different SMBH mass distributions (compare Fig. [7] with 
[5]). Consequently it is better not to go on easy conclusions 
regarding the activity and the evolution of both AGN and 
normal host galaxies, until we have understood what the 
best scaling law able to infer correctly the SMBH masses 
is. 
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A Log(R e \cr) A Log(M.|tr) 

Fig. 6. Correlations b etween the residuals in the M, — a and R c — a (left), and M* — a (right) relations at each u 
(Hop kins et a l. 2007b), for the full sample of 57 galaxies (a-b), and for a more consistent sample (see text) of 52 galaxies 
(b-c). The corresponding intrinsic scatter and Pearson linear coefficient are reported in the upper-left corner of each 
plot. The symbols are the same as in Figure [TJ Galaxies inside the black circles are pseudo-bulges. 



6. Conclusion 

We analyzed different scaling laws for a consistent sample 
for med by N = 52 galaxies, which have been catalogued 
bv lSani et ail (|2011[ ) on the base of Spitzer /IRAC 3.6 /xm 
observations of local Universe. The sample is formed by 
both early-type and late-type galaxies. Actually, from the 
original sample of 57 galaxies, w e removed 5 disc galax- 
ies identified bv lSani et al.l (|201lh as hosting pseudobulges 
and that are non consistent with the correlations for classi- 
cal bulges. For the galaxy masses, we considered both the 
dynamical mass and the stellar mass. The results of our 
analysis have been reported in Tabled and Figures 1-4. 

The main emerging result is that the relation between 
the mass of SMBHs and the kinetic energy of random mo- 
tions of the host local galaxies appears to be a robust corre- 



lation, which could provide the right passkey to understand 
the nature and evolution of the numerous observed corre- 
lations between SMBHs and host spheroid properties. This 
is in agree ment with our pre v ious studies performed o n the 
sam ples o f[Graham| (|2008aD . iGiiltekin et ail ((2009a). and 



|2009[) . Even if the values of the Xi (see Tabled]) indi- 



cate that the M, — Mqct 2 works better than the others, by 
considering the intrinsic scatter of the various relations and 
in particular their errors, we cannot conclusively determine 
the best relation, because all the examined laws appear 
on the same level. The comprehensive analysis of residu- 
als discussed in flUdq es not confirm the result claimed by 
Hop kins et al.l i|2()0 7l>; , according to which a black hole FP 
should be preferred over a simple one-one relation, even if 
it cannot be definitively ruled out. 
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Fig. 7. The distribution of SMBH masses for both normal (white) and AGN (colored) host galaxies as a function of 
morphology (A[log 10 M,j = 0.05 bin). The data of the galaxies have been extracted from the Sloan Digital Sky Survey 
in redshift interval 0.02 < z < 0.05. The masses have been inferred via the M, — a relation using t he slope and the 
normalization taken from Tabled In order to compare our results with that of iSchawinski et all (|2010[) . we plotted only 
objects where the measured velocity dispersion is greater than 40 km s _1 . 



Since it has now been tested on four different samples 
of galaxies, independently catalogued by different authors, 
the goodness of the M, — Mqct 2 relation as a predictor 
of the SMBH mass in the center of galaxies is enough ro- 
bust. Again, this relation is the only one that currently has 
a quite clear physical explanation. In this perspective, as 
we discussed in JJSJ in order to obtain correct estimates of 
SMBH masses, the M. — Mqct 2 relation should be prefer- 
ably used instead of the other popular scaling laws. As a 
matter of fact, the SMBH mass distribution of the early- 
type AGN galaxies, inferred by a much more physically mo- 
tivated relation, clearly shows the radio-quiest/radio-loud 
dichotomy. The same result is not achieved if we use the 
usual Mm — a relation. 
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